RAL-TR-98-042 



hep-ph/ 9805111 



Secondary Heavy Quark Pair Production 
in e+e~ Annihilation 

D. J. Miller 

and 

Michael H. Seymour 

Rutherford Appleton Laboratory, Chilton, 
Didcot, Oxfordshire. 0X11 OQX. England. 



Abstract 

The multiplicity of heavy quarks from gluon splitting in e"'"e~ annihilation 
has now been theoretically calculated and experimentally measured at LEP. 
However, the experimental measurement requires theoretical input for the shape 
of the multiplicity with respect to an event shape. In this paper we calculate 
the multiplicity of heavy quarks from gluon splitting in e+e" annihilation as 
a function of the heavy jet mass of the event, to next-to-leading logarithmic 
accuracy. We compare our result with Monte Carlo predictions. 
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Introduction 



Heavy quark production in e"'"e~ annihilation can come from two sources: from the hard 
interaction itself, e"'"e~ QQ, and from the splitting of perturbatively produced gluons, 
e+e~ — > qqg — > qqQQ. We call the latter secondary heavy quarks. Their total rate is an 
infrared-safe quantity, so can be calculated as an order-by-order perturbative expansion in 
as, starting at (9(ag). The leading order term in this expansion was calculated in Refs. [^0]. 
At higher orders in as, large logarithms arise, ag log^"~^(s/mQ), potentially spoiling the 
convergence of the perturbative series at high energies, s ^ uiq. In Ref. P] these (lead- 
ing) logarithms were summed to all orders in as. The next-to- leading logarithmic terms, 
ag log^"~^(s/mQ), were resummed in Ref. ^ yielding a result that is uniformly reliable for 
all s. However, only the total multiplicity was calculated, retaining no dependence on the jet 
kinematics. In this paper we perform a more differential calculation, allowing the multiplicity 
to be calculated for various event shapes. 

Several experimental measurements of the secondary heavy quark production rate have 
now been made. In Ref. it was extracted for charm quarks from a measurement of 
the D* fragmentation function, and found to be more than a factor of two above the ex- 
pectation of Ref. although with large systematic errors coming from uncertainty in the 
fragmentation function of primary charm quarks. References ^-0] made less model depen- 
dent measurements by selecting hard three-jet events, which enhances the fraction of heavy 
quarks produced by the gluon splitting mechanism. In general the measurements have been 
above the prediction of Ref. [T[] , although within the range allowed by variations in as and 
the quark mass. 

For the parameters 

as = 0.118, (1) 

TJic = 1.2 GeV, (2) 

mb = 5.0 GeV, (3) 

we obtain results for the fractions of Z° decays that contain a secondary charm or bottom 
quark pair, of 

/c = 2.007%, /b = 0.175%. (4) 

Notice that these are somewhat higher than those of Ref. [Q], principally because we have 
used a different value for asQ. These should be compared with the LEP values^ of 

/c = (2.44 ± 0.43)%, /b = (0.22 ± 0.13)%. (5) 

It should be noted that the summation of the large logarithms is essential for this agreement. 
For example, the leading order result alone for fc is only 0.878 for the same parameters. 

However, these measurements are still somewhat model dependent, since they rely on 
Monte Carlo event generators to correct from the rate seen within the selected three- 



jet region, to the total rate. Since in this paper we calculate this dependence explicitly, it 

*For our parameters, the calculation of Ref. gives fc = 2.200% and /b = 0.207%. The remaining 
difference is due to next-to-next-to-leading logarithms. 

^For the bottom quark there has been only one measurement |^ , while for the charm quark, we have 
averaged the results of Refs. and 0, assuming that the systematic errors are uncorrelated. 
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will become possible to directly compare experiment with theory in the measured region, 
without the need to extrapolate to the unmeasured region. 

We begin by outlining the calculation of the resummed multiplicity, retaining the de- 
pendence on the jet kinematics, and thereby allowing the calculation of any event shape 
variable to be performed numerically. We then present numerical results for a few quantities 
of interest. Finally, we make some comparisons with the Monte Carlo event generators that 
have been used previously, and draw some conclusions. 



Calculation of the Resummed Multiplicity 

Fixed Order 

In order to separate the processes of primary and secondary heavy quark production, it is 
necessary that the interference between them be zero, or at least very small. Fortunately 
this is the case[0. For non-identical quarks coupled to a vector current, the interference term 
will vanish by Furry's theorem if we assume that the charges of the quarks are not measured 
(i.e. we do not distinguish quarks from antiquarks). For an axial current, cancellations will 
occur between up- and down-type quarks, leaving only the case where the "light" quark is a 
bottom quark. This will only provide an effect of the order of 0.2% of the secondary heavy 
quark rate[jll|. The contribution for identical quarks is only slightly larger. In this case 



Furry's theorem no longer applies but necessarily the quarks will be of the same flavour - 
either bottom or charm, depending on the case in question, and can be neglected. 

Then, the leading order differential cross section for the production of secondary heavy 
quarks, 7* —>■ qqQQ, is easily calculated. In Ref. |Q, this expression is integrated over the 
heavy quark momenta, yielding the simple result given in Eq. (2) of Ref. In our case, 
since we must retain the jet kinematics, we cannot do this integration since we would not 
correctly account for events in which the heavy quarks fall into different hemispheres or 
different jets. Instead we must retain the full unintegrated amplitude, which does not have a 
compact form. This expression is then integrated over the appropriate phase space by Monte 
Carlo methods in which we can apply any event shape definition desired. 

We choose to calculate the resummed multiplicity of secondary heavy quarks as a function 
of the heavy jet mass. Here we use the thrust definition of heavy jet mass, i.e. we separate 
each event into two hemispheres via the thrust axis and take the heavy jet mass to be the 
larger of the two hemisphere masses. Our calculation is therefore easily converted to give 
the differential distribution with respect to any thrust-like event shape. 



Resummed Calculation 



There are two requirements to which we must conform in our calculation of the logarithmic 
contribution. Firstly we must retain the exact kinematics of the qqg production in order to 
be able to accurately obtain the heavy jet mass for each event. We need not worry about 
the exact kinematics of the heavy quarks, since the large logarithms arise from the parts 
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of phase space where they become colhnear and only the kinematics of the gluon need be 
considered (although the exact kinematics of the heavy quarks are, of course, included in the 
fixed order contribution). Secondly, we must include all soft gluon emission from the light 
quarks and virtual gluon. These emissions contribute to the heavy jet mass (by making the 
light quarks massive) through large logarithms, which must be summed to all orders using 
the coherent branching formalism. 

Bearing these considerations in mind, we take the differential multiplicity to be, 

^QQ_ (^Mff, Q ; Qg) = J dXi dX2 dki dkg ^2 ^Q2 ^9(^1? ^lmax)fq{^2-} ^2max) 

nf{kl kl; Ql)e{kl - Ql)5{Ml - /.(xi, x^, kl kl k])). (6) 

The above notation requires some clarification. As usual, x\ and X2 are the energy fractions 
of the light quark and antiquark respectively, and fci, k2 and kg are the four-momenta of the 
quarks and gluon. The maximum value of kf^ % = 1,2, as constrained by the phase space 
limits is kf^^^^. Also, kj_ is the transverse momentum (squared) of the virtual gluon, given 

by, 

fci = (1 - xi + ei - e2)(l - X2 + 62 - ei)g^ (7) 
where is the centre-of-mass energy squared and 
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mf , , 

^^=gJ' ^ = 1'2 (8) 

are the rescaled (primary) quark masses (squared), which result from the soft gluon emission. 
Here, h{xi,X2, kf, k^, kg) is the heavy jet mass as a function of the exact kinematics of the 
quark, antiquark and gluon. 

Summation to all orders of the leading and next-to-leading logarithms is included in the 
functions fg and n^^. The function fg is the quark jet mass distribution which has been 
calculated to next-to-leading logarithmic accuracy in Ref. [|12|. More exphcitly, fq^k"^, Q^) dk"^ 
is the probability that a quark created at a scale gives rise to a jet with mass squared 
between k"^ and k^ + dk'^. This function includes all soft gluon emission from the light quarks 
and sums to all orders leading and next-to- leading logarithms of kf^^^/kf, i = 1, 2. 

The function n^^ is the gluon jet mass distribution weighted by the heavy quark pair 
multiplicity, and has not been calculated elsewhere. In other words, n'^'^{kg, kj_; Ql) dkg is 
the number of heavy quark pairs within a gluon jet that was formed at a scale k"^ and has 
a mass between kg and kg + dk"^. In calculating this quantity we must be sure to include 
leading and next-to-leading logarithms of both k'^/kg and kj_/Ql. 

Notice that the resolution scale at which the heavy quarks are resolved is, 

go = 2m^ = mQe^/^ (9) 

Here, rriQ is the heavy quark effective mass. It has been shown in Ref. that by using 
this effective mass we can neglect the heavy quark mass in the decay of the gluon while 
maintaining the correct leading and next-to-leading logarithms]^ 

■l-This follows from the comparison of Eqs. (16) and (17) of Ref. ||l|. 
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The Multiplicity Weighted Mass Distribution 

As in the case of the jet mass distribution, it is more convenient to calculate the integrated 
distribution, 







(10) 



Physically this is the number of QQ pairs resolved in gluon jets of mass squared less than k"^. 
It can be derived from N^{k'^,Q'^;Ql), the multiplicity of gluons within gluon jets of mass 
squared less than k"^, by integrating over the kernel for the splitting g QQ, i.e. Pqg. There- 
fore, we have, 



dz^P,,iz)N^^ik\Q';q') 

m{k\Q'-q'). 



1 r^^ dq^ asiq) ^ 



has been derived in Ref. [0, and is given by. 



where. 



N^^{k\ Q'; Ql) = F,{k\ Q') {Ml{k'; Ql) + {P{k'; Ql) - P{k'/Q^; Ql)) } (12) 



-iC (2^/{zo,Zk)-zl/zl + l 



and J^-j is the usual multiplicity of gluons within a gluon [114 



Af„'ik';Ql) = M+{z,,Zk)-CM{z,,zu). 



(13) 



(14) 



Here, we have used the variables introduced by Catani et al.||Tj], where. 



2,2'kCa 



The functions, A/"^ and Af, are defined in terms of Bessel functions by. 



J\f^{zo,Zk) 
J^{zo, Zk) 



Zk 
B 



B 



[lB+liZk)KB{zo) + KB+l{Zk)lB{zo)] 



[lB{Zk)KB{zo) - KB{Zk)lB{zo)] 



and the parameters B and C are given by, 

1 /II 2Nf ACpNp 

^-bKY^^^S scT 

where b is the first coefficient of the /3-function, 



b = —Ca - -Nf. 



3 b Ca' 



(15) 



(16) 



(17) 



(18) 



The gluon jet mass fraction, Fg{k'^, Q^), is the probabihty that a jet formed at scale will 
have a mass of less than k"^. This is the integrated version of fg{k'^,Q'^) and is given, in 



Ref. m, by 



F.ik^Q') 



expjcAlogf^j /i 



' as{Q) 
Arc 



where, 



/i(A) 
/2(A) 



2 

bX 
2 



^log 



(1 - 2A) log 
1 



1 - 2A 



2(1 - A) log 



1- A 



1 - A 



(19) 

(20) 
(21) 



Notice that the result of Ref. |T2| includes leading and next-to-leading logarithms within 
logFg, whereas the expression above is accurate only to leading and next-to- leading loga- 



rithms within Fg itself. 



The integrals required for the evaluation Eq. (|TT]) have been derived in Ref. [0] and are, 
to next-to-leading logarithmic accuracy. 



dz 



zo z 

dz 



— Af+{z,Zk) = Af,{zo,Zk)-^{l-B{l-^]] ^+{zo,Zk) + ^, (22) 



J\f{z, Zk) 



■ Af' {Zq, Zk) - \ 



'20 Z ZqZi 

This gives, 

Nfik', Q'; Ql) = F,{k\ Q') [ Mf{k', Ql) + {l<^%k', Ql) - I^^ik^Q^; Ql)) } 



(23) 



(24) 



with, 



3b C A 



(25-2-C)-i(C + 2)log 









c ' 




ZoZk 




1 Zk 





U+{zQ,Zk) 



(25) 



and. 



J^f{k': Ql) = irA ^lizo, z,)+ -{B-1 



3b 



ZoZk 



-(25 + C)j Af+{zo, Zk) + _(C + 2)| . 

(26) 

In the above and following A/i is defined as A/" but with B replaced by 5 — 1. 

Finally n'^'^ is obtained by differentiation with respect to the jet mass, yielding the result, 

nfik^Q';Ql) = y,Fg{k\Q') [M'^^^ik'-^Ql) + {l'^^ {k\ Ql) - 2 {k^ Ql))} 



+ ^,{\ogFg{k\Q'))Nf{k\Q'-Ql), 



(27) 



with, 
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60 \Zq ZqZi I 



-fi(i-4/.o=)). 



Af{zo, Zk) 

{2i 



and, 



A = ^Mog(QA). (30) 

Notice that the matching to tree level is particularly simple. We only require the expan- 
sion of n'^^ to order as? which is, 

nf{k\Q'-Ql) = ^^^ + 0{al). (31) 



where, 



Calculation of the Background 

The background to secondary heavy quark production in e"'"e~ annihilation, i.e. primary 
heavy quark production, can be estimated by standard three-jet production, since the mass 
effects will be small. For the analysis of the heavy jet mass, the fixed order contribution to 
this background is given by, 

1^ = 1 C,^j-41ogr-3 + 6rlogr-4-^log(l-r) 

ao dMfj Mfj 271 y 1 — r 

-F— ^ log (1 - r) - 6r log (1 - 2r) - log (1 - 2r) + gr^l (32) 

1 — r 1 — r J 

where cx*^^-* is the 0{a^) contribution to the cross section, ctq is the Born cross section, 
r = Mfj/Q"^ and is the heavy jet mass. 

Of course, the large leading and next-to-leading logarithms must also be included and 
are given by, 

—^-Jm^ = ^f^i^H, Q')F,{Ml, Q\ (33) 
Expanding this to 0{a^) gives the first two terms of the fixed order piece. 

Therefore it is clear that matching with the fixed order will result in the full answer, 

= 2/,(M^,g^)F,(Mi,g^) + -^C^^{6rlogr-4-^log(l-r) 
(Jo aMfj Mfj 27r [ 1 — r 

^ log(l -r) -6rlog(l -2r) ^ log (1 - 2r) + Or^j . (35) 



1-r °^ ' . 
This can be seen plotted in Fig. |l[ 
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Numerical Results 

For all the distributions we show, we concentrate on their shape, normalised to the number 
of secondary heavy quarks, rather than on the total rate. We use the as and quark mass 
values quoted earlier. 

We present the heavy jet mass distribution for a/s = mz in Fig. [l|. This is closely related 
to the jet mass difference, Mh — Ml, which was the event shape used to fit fc in Ref. 0. We 
see that the heavy jet mass provides a good discriminator of events with secondary heavy 
quarks from the three-jet background. 




Figure 1: The multiplicity of heavy quark pairs in decays as a function of the heavy jet 
mass, normalised to the number of secondary heavy quarks. The shapes for b quark (solid 
curve), and c quark (dashed) pairs are compared to the three-jet background (dotted). It is 
clear that secondary heavy quark production can be distinguished from the background by 
the shape of the heavy jet mass distribution. 

Of course, these shapes are dependent on the values chosen for the parameters, Aqcd 
and niQ. The effect of varying these parameters is seen in Fig. 0. Also shown is the contri- 
bution from the fixed order term alone, demonstrating the importance of resumming large 
logarithms. 



Event Generators 



Monte Carlo event generators predict quite a wide range for the rate of secondary heavy 
quark production. While JETSET ^ and HERWIG [Q, which are both based on the parton 
shower formalism, are in quite good agreement, ARIADNE [1^, based on the dipole cascade 
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Figure 2: The multiplicity of bottom quark pairs as a function of the heavy jet mass, 
normalised to the number of Z° decays. The dashed curve shows the fixed order result, 
whereas the solid curve includes the resummation of the large logarithms. Also shown are 
the results of varying the quark mass by 5% (dotted) and Aqcd by a factor of two (i.e. as 
by 10%) (dash-dotted). 

model, lies well above them. At present the data for the total rate lies between the two, in 
agreement with both, although in somewhat better agreement with ARIADNE. 

In Ref. [H] it was argued that the discrepancy between the models is actually due to a 
specific problem with ARIADNE - the fact that it allows very low transverse momentum 
gluons to be very massive. If that is the case, it should show up in the distributions calculated 
in the previous section. Following the suggestion of Ref. , later versions of ARIADNE have 
had an option^ to veto gluon splitting with rrig > kj_g, which should fix this problem. 

The parton shower and dipole cascade models are only formally accurate to leading 
logarithm and do not include the exact matrix elements for qqQQ production. Therefore 
our calculation is more accurate than them and can be used to check them. 

In Fig. I we compare our results with the predictions from the event generators for 
a/s = mz- With the exception of the gluon splitting option in ARIADNE, we keep all model 
parameters at their default values. We see that HERWIG and JETSET give similar predic- 
tions for the distribution as well as for the rate and that the unmodified ARIADNE peaks 
at somewhat lower heavy jet mass than them. Adding the new modification, ARIADNE's 
distribution is more like the other models', but still somewhat different, particularly at low 
jet masses. Our results lie between ARIADNE and the other models. 

Increasing the centre-of-mass energy, the relative importance of the fixed order term 
§This option is switched on by setting MSTA(28)=1. 
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Figure 3: The multiplicity of bottom quark pairs as a function of the heavy jet mass in Z° 
decays from our calculation (solid) and from various Monte Carlo models. 

is reduced and one gets a cleaner probe of the parton evolution. In Fig. ^ we show the 
comparison again at higher energy. The modified version of ARIADNE is in even better 




Figure 4: As Fig. | but at y/s = 500 GeV. 



agreement with the other two models, while the unmodified version is in good agreement 
with our calculation. We therefore see no evidence to support the claim of Ref. |jl|] that there 
is a problem with ARIADNE. 
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Summary 



We have calculated the multiplicity of heavy quarks from gluon splitting in e"'"e~ annihilation, 
as a function of the heavy jet mass. Our result is exact to leading order in ag, and sums 
leading and next-to-leading logarithms to all orders in as- We find that the fractions of Z° 
decays that contain a secondary charm or bottom quark pair respectively, are 

/c = 2.007%, /b = 0.175%. (36) 

The shape of our result is similar to that predicted by Monte Carlo event generators at the 
Z°, lying between the different models, but in better agreement with ARIADNE at higher 
energy. 
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